Relaxation to equilibrium of the expectation values in macroscopic quantum systems 
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A quantum mechanical explanation of the relaxation to equilibrium is shown for macroscopic 
systems for nonintegrable cases and numerically verified. The macroscopic system is initially in 
an equilibrium state, subsequently externally perturbed during a finite time, and then isolated for 
a sufficiently long time. We show a quantitative explanation that the initial microcanonical state 
typically reaches to a state whose expectation values are well-approximated by the average over 
another microcanonical ensemble. 
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I. INTRODUCTION 

Q ' It is our experience that isolated thermodynamic systems would eventually relax to a stationary state called 

O \ equilibrium [l|, y]. In thermodynamics, this statement is regarded as requirements for the relevant choice of the 
model and its macroscopic descriptions. Recently, it is shown that equilibrium states are very popular in the ex- 
I I tremely large Hilbert space of state vectors for macroscopic quantum systems [3-[l2|. The expectation values with 
d • respect to a randomly sampled pure state are typically well described by equilibrium ensembles [3, |^, [l2|. And it is 
ti ■ expected that initial nonequilibrium state approaches to typical state of equilibrium. Then it is reasonable to explore 
a dynamical explanation for the macroscopically observed behavior of the relaxation to equilibrium [J, |5|, |8|. Indeed 
Cd ' the coarsgrained macroscopic quantities typically yield a decomposition of the total Hilbert space with a subspace of 
overwhelming dimension which represents equilibrium[4!, '5!]. Based on this property, arbitrary state approaches and 
stays near equilibrium[4-6], which is a modern version of quantum ergodic theorem[13j. By assuming the thermo- 
' ^ ' dynamically normal fluctuation of the quantities, the approach to equilibrium can be derived[ll|. Also the reduced 
density matrix of huge isolated system approaches to time averaged state and stays near equilibrium for almost all 
the time|8j]. Description of the relaxation is also important in the context of nonequilibrium theorems |lj,|l5|, which 
relies on the use of thermodynamic ensembles. As for the isolated system, the relaxation is not guaranteed for the 
total density matrix itself. Indeed the total density matrix does not actually relax to a thermal state in the course 
of unitary time evolution. Linear phonons in the initial microcanonical state which are externally perturbed and 
excites only few modes will not approach to a microcanonical state. For classical systems, there exist integrable 
systems which do not reach an equilibrium state even after a long waiting time due to the conserved quantities. For 
isolated quantum systems, there are always conservative quantities such as projection operators on a specific energy 
eigenstatesp, [ifl, H; ll3| • And it has been investigated which conservative quantities would affect the thermalization 
r — I . property [l6|. The thermalization of a few body quantities has been established for some nonintegrable systems whose 
^^ ' energy spectrum do not tend to degenerate [13, 13 • There are also reports that the reduced equilibrium ensemble 
average works for spin chains both for integrable and nonintegrable parameters'19], which is consistent with the typi- 
cality of thermal statc[3] but in marked contrast to jl7|,il8|| at least for the small system sizes. It is possible that there 
are formulations of relaxation properties which are not very sensitive to the nonintegrability. Since we assume the 
initial microcanonical state, thermalization for the expectation values of the conservative quantities is trivial for au- 
tonomous systems as explained later on. Also the evaluation of off-diagonal elements for the nonintegrable case shows 
quantitatively good agreement with the numerical values, which is a sufficient condition to guarantee the relaxation 
property. 

We insist that from the uniformly random sampling assumption which holds for nonintegrable cases and the initial 
equilibrium state, the expectation values of quantities which polynomially depend on the system size typically show 
relaxation to the equilibrium values. Especially, it is pointed out that the off-diagonal elements of the density matrix 
does not contribute to the expectation values. In this article, we explain how the statement is made quantitative 
for the time-dependent unitary time evolution. And the explanation is numerically verified for a spin chain with a 
magnetic field in the nonintegrable regime. This paper is organized as follows. In Sec. II, we present our model with 
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the precise external forcing procedure. In Sec. Ill, the relaxation of the expectation values are explained. Sec. IV is 
devoted to a summary. 

II. MACROSCOPIC SYSTEMS 

Let us consider a macroscopically large but finite initially isolated system. After initial time t = 0, an external 
forcing acts on the system, and the total Hamiltonian H{t) depends on time. Until t = 0, the density matrix describing 
the initial state is assumed to be microcanonical specified by an energy E 

pm-^^siE-nm, (1) 

where J7(0) is the density of the states. In the course of time evolution, the external work is done on the system 
through the time dependence of H{t). 

The deterministic external forcing acts during the time interval < i < Tq, and switched off for t > Tq. It is 
expected that after a sufficiently long waiting time, i.e. at t = T ^ Tq, the density matrix reaches a state which 
yields approximately the same expectation values as the microcanonical ensemble p{T) = ^^}j'\ 5{E + AE' — H{T)) 
with an energy shift Ai?. 

III. RELAXATION OF THE EXPECTATION VALUES 

It is our statement that under the uniformly random sampling assumption valid for the nonintegrable cases, the 
microcanonical ensemble well approximates actual expectation value of quantity A whose maximum eigenvalue poly- 
nomially depends on the system size 

TtUp{Q)U+A = Trp( T)A. (2) 

This evaluation shows that the relaxation to equilibrium is explained as the property of the expectation values instead 
of the density matrix. Note that the initial microcanonical state guarantees that Eq.® holds for the conservative 
quantities with p{T) = p(0) for autonomous systems irrespective of the integrability. An important class of observables 
is extensive quantities which are sums of local quantities such as total magnetization in a spin chain. Such quantities 
mutually almost commute, since the local quantities at spatially different points commute. 
The actual final state is reached by a unitary time evolution 

Up{Q)U+ 

= 2_^^n,m\En){Em\, (3) 

where we presented the eigenstate of H(T) as H{T)\En) = En\En), Cn,m = 

e-T^(E^-E^)('^-'ro)(^En\U{To)p{0)U{To)+\E^), and U{To) = T{e- fo" jiHit)dtj^ After t = Tq, the total system 
is isolated, and evolves by the Hamiltonian H{T). On the other hand, for the nonintegrable systems the off-diagonal 
elements of the physical quantity A is negligible compared to the diagonal elements: 

i) The off-diagonal matrix element {En\A\Ejn) of the nonintegrable system is typically negligible in the macroscopic 
limit. Physically, this would be reasonable since it means that transition amplitudes between macroscopically 
different states due to the perturbation A are extremely small[20|. We shall give a quantitative explanation of 
this statement. The analysis is based on the high-dimensionality of the Hilbert space as in Refs.|3l.[7l.[9l.[l2l|. Let 
us diagonalize the quantity A as 

A = Y,An\A^){Anl (4) 

n 

and define its square root 

i3 = vyi;i^„)(A„|. (5) 



Here the spectrum {An} is assumed to be nonnegative. Without loss of generahty, the operators bounded below 
such as energy and number of particles can be made into this form. We also define state vectors {|$„) — B\En)} 
so that the matrix element is expressed as the inner product 

{E.r,\A\Ern) = ($„|$,„). (6) 

The states {!'&«)} are chosen from the extremely large Hilbert space H with various directions. Thus we assume 
that the sequence of the normalized vectors { , , ...} is regarded as a uniformly random samplinq 

from divaH — 1 dimensional unit sphere as we will confirm for the case of nonintegrable spin chain. It is then 

straightforward to show that the mean square of the inner product is smaller than -r^-'W [3] , which we will show 
later. Here || A|| is the maximum of the eigenvalues of A. 

Let us derive the inequality for the inner product 

(l('i'n|$™)P)< j^, (7) 

where the bracket shows the average with respect to the uniform random sampling of , ' "^ from the unit 
sphere. Uniform random sampling is expressed by the vector representation of . ' "' in an orthogonal com- 
plete basis as (i„ = (cos 9i , sin 9i cos 6*2 , sin 9i sin 02 cos 63, ..., sin 9i sin 62 ■ ■ ■ sin 6d- 1 ) with d — dimH . Hereafter 
we use the abbreviated notation for the dimension d. The angles 9i are uniform random variables on the real 
axes. We are interested in particular cases where \\A\\ polynomially depends on the system size. On the other 
hand, the dimension dimH grows exponentially as the system size increases. Then inequality ([7]) shows that 
the off-diagonal elements Eq.® is extremely small for the macroscopic system. We numerically show that 
Eq. ([7]) actually holds and hence random sampling assumption is reasonable for a quantum spin chain. The spin 
chain is regarded as a nice example, since it provides observables, which is a sum of the local quantities. The 
Hamiltonian of the spin chain is chosen as 

N N N 

where J stands for the exchange interaction between neighboring sites, and a constant magnetic field B = 
(a, 0, 7) is applied, cr* is the i = x,y, z component of the Pauli matrix at the site j. 7 controlcs the nonintegra- 
bility. As a macroscopic quantity, we consider the square of the x component of the total magnetization 

N 

Other choices of A are possible, but Eq.® guarantees a well-defined square root B = J2j=i '^j- I^O'^ example, 
if we choose A as x component of the total magnetization 01 A = {J2T=i '^'jY "* — ^j ^^ still observe behavior 
resembles to Fig.l. The system sizes 6 < A^ < 11 are explored. Note that the dimension of the total Hilbert 
space is 2^ and can be large for relatively small N . By diagonalizing the Hamiltonian matrix H under the 
periodic boundary condition, numerical eigenvectors \En) are obtained, and the absolute values of all the matrix 
elements |(£'m|^|i?n)| are shown in Fig.l. The data are shown only for parameters J = 1, a = 1, and 7 = 0.5, 
however, we have confirmed similar behavior for various values of 7 > 0.1. To verify that Eas. (j7ll0|) are satisfied 
quantitatively, we also show the matrix elements |(£'iooo|^|-E'n)| with a fixed value of ?7i = 1000 for N — 11. 
The statistical mean value of the off-diagonal elements agrees with the theoretical estimation -^ = -ILJ- as 
shown in Fig. 2. This estimation is considered as a generic property of macroscopically different states for the 
nonintegrable systems, since it is expressed only by the typical value of A and the dimension d. Thus for large 
d = 2^, the directions of |$„) are orthogonal to each other, and are considered to be distributed on the d— 1 
dimensional unit sphere without any bias to a specific direction. 
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FIG. 1: The absolute values of matrix elements |(_Em| j4|i?„)| of ^ = {^ ^ 



The average of the off-diagonal elements is the 



same order as the average of the diagonal elements divided by \/d. The system size is A'^ ; 
similar results for 6 < A'^ < 11. 



We have confirmed qualitatively 
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FIG. 2: The absolute values of matrix elements |{_Eiooo|^|S„) of A 
The red line shows average of the off-diagonal elements {|{-Biooo|^|£'n)|) 
diagonal elements divided by ^2^, 0.138. These two lines are mutually very close. 



,^j^-jy . y and z components give similar results. 
- 0.173, and the black line gives the average of the 
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3: The absolute values of matrix elements |(i?iooo|^|£'n)| oi A = (^D^^i ''"J )^- V s-^d z components give similar results. 
The data are shown for parameters J = 1, a = 1, and 7 = 0, where the Hamiltonian is integrable. The mean of off-diagonal 
elements is 0.03719, and the theoretical value is K-Eioool^iooo) ^ 0.1068. 



is immediately derived, since the expectation values of each component are the same and their sum should be 
unity. Since \En) and l-Bm) are normalized, 

($„|$n)(*m|$m) 

< Uf. (11) 

The inequality (O derives from Eqs. (|10|lip . For comparison, the matrix elements |(_Eiooo|^|^n)| at A^ = 11 
for the integrable case 7 = are shown in Fig. 3. As in Fig. 2, A = {J2T=i ^]Y is explored. For m < N and 
i = y,z, \{EiQQQ\A\En) behaves as in the nonintegrable cases. But for i — x, the smallness of the off-diagonal 
elements is affected. Aside from the dominant contribution from the diagonal element, there are subdominant 
contributions as shown in Fig. 3. The same tendency takes place also for m < N. This shows that the uniformly 
random sampling assumption is affected at least for the relatively small system size. 

From Eq.©, the coefhcient c„_m of Eq.® is safely replaced by {En\U p{0)U~^\Ejn)Sn^m in the evaluation of 
TrUp{0)U+A = Y.^^jE^\Up{6)U+\En){Eri\A\E^). The off-diagonal elements of Up{0)'u+ does not contribute 
to the expectation value of the quantity A, since |(i?m|^|£^n)| {m ^ n) is neglected.. More quantitatively, the 
error caused by neglect of the off-diagonal elements is indeed negligible by random phase approximation, since 
there are (P terms of O(^) with various phases. 

ii) The diagonal elements of the density matrix Up{0)U^ is wcU-approximatcd by those of a microcanonical state 
with respect to H{T). It is important to note that the initial microcanonical state specifies an energy scale E. 
Indeed the diagonal elements are expressed as 

{En\Up{0)U+\En) 

{E„\S{UH{0)U+-E)\Er,) 



n{E{0)) 



= n(^)T.^(^^rn-E)\{E^\E^)\^, (12) 

where we introduced the normalized eigenstates of UH{0)U^ as UH{0)U^\En) — E,i\En), i.e. U\En{0)) — \En) 
with an eigenstate oi H{0), |ii^„(0)). The set of states {|£'„)} and {|£'„)} are related by a unitary transformation 
J2n \En){En\- In the presence of perturbation, |(£'„|i?m)P — \{En\U\Em{0))\-^ would be non negligible only when 
the conservation of the energy is well-satisfied after the long waiting time, i.e. En = Em. + A_E with the energy 
change AE from the initial to final times caused by external perturbation during < t < Tq. 

We show a quantitative estimation of AE. The matrix element (i?„|_E,„) is evaluated as 6nm + ^-^ — l (1 ~ 
i^Tim) up to the first order of the perturbation H{T) — C/iJ(0)C/+, where the factor (_E,„|i7(r)|_E„) would be of 



order -J - — ^ , '" as in Eq.©. Then it is immediately shown that Em contributes to {En\Em) only when Em and 
En are sufficiently near \Em — En\ < -y^, otherwise the ratio ^ _1 is negligible. Such \Em) is expanded 

as \Em) = X^n ^n,nil^n)j whcrc the coefficient dn,m is non negligible for £"„ sufficiently near Em- The energy 
change AE is determined as 

AE = {Em\{H{T) - UH{0)U+)\Em) 

= 2^Mfem| Ek — Em 
k 

= En-E, (13) 

which is almost independent of the suffix n. Here we evaluated as J2k \dk,m\'^Ek = Em = En and Em — E 
from the Dirac delta in EQ. (fT2|) . The important property used here is a continuity of the mapping, i.e. when 

\En — Em\ is small enough compared to y ^, \En ~ Em\ is also sufficiently small. 

Therefore the third line of Eq. ([T2l) has a sharp peak at En = E + AE as a function of En and is proportional 



to the function S{En — E — AE) as 



{En\Up{0)U+\En) 
^^SiE^-E-AE) 

iEn\-^S{H{T) - E - AE)\E^), 



(14) 



where the density of the states at i = T is determined uniquely from the normaUzation, and H(T)\En) = £'„|£'„) 
is used. The diagonal elements is thus given by the microcanonical ensemble. 

Then as far as the expectation value is concerned, the state Up{0)U^ should be well-described by the microcanonical 
ensemble p{T). 

IV. SUMMARY 

In conclusion, by assuming that the initial state is prepared as a microcanonical ensemble and the uniformly 
random sampling assumption holds as in the nonintegrable systems, the actual state at time t — T can be replaced 
by another microcanonical ensemble in the evaluation of quantities which polynomially depend on the system size. 
The derivation is based on the high dimensionality of the Hilbert space and consequent uniformly random sampling 
assumption valid for the nonintegrable systems, and restrictions for the system size dependence of observables as 
well as the perturbative treatment of unitary transformation of the energy eigenstates. The validity of the uniformly 
random sampling assumption for . is numerically verified for a nonintegrable spin chain. It is also remarked 

that the initial state can be out of equilibrium, i.e. the initial density matrix is p{0) — UojT77yrS(E — H{0))Uq with a 

unitary transformation Uq, since this can be regarded as the state evolved from an actual initial state Uq p(0)Uo = 
^S{E-H{0)). 
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